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ABSTRACT 

We study the effect of short term variations of the evolution of AM CVn systems on 
their gravitational wave emissions and in particular LISA observations. We model the 
systems according to their equilibrium mass-transfer evolution as driven by gravita- 
tional wave emission and tidal interaction, and determine their reaction to a sudden 
perturbation of the system. This is inspired by the suggestion to explain the orbital pe- 
riod evolution of the ultra-compact binary systems V407 Vul and RX-J0806+1527 by 
non-equilibrium mass transfer. The characteristics of the emitted gravitational wave 
signal are deduced from a Taylor expansion of a Newtonian quadrupolar emission 
model, and the changes in signal structure as visible to the LISA mission arc deter- 
mined. We show that short term variations can significantly change the higher order 
terms in the expansion, and thus lead to spurious (non) detection of frequency deriva- 
tives. This may hamper the estimation of the parameters of the system, in particular 
their masses and distances. However, we find that overall detection is still secured as 
signals still can be described by general templates. We conclude that a better mod- 
elling of the effects of short term variations is needed to prepare the community for 
astrophysical evaluations of real gravitational wave data of AM CVn systems. 
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1 INTRODUCTION 

The Laser Interferometer Space Antenna (LISA) is an 
ES A/NASA space-based gravitational wave (GW) laser in- 
terferometer designed to observe a wide range of sources in 
the frequency range ~ 10"^ - 1 Hz (|Bendeill998h . At the mo- 
ment LISA is in its design stage with an open launch window 
of 2018+. One of the important classes of objects are Galac- 
tic compact object binary systems, in particular double 
white dwarf binaries (|Evans et al]|l987l : iLipunov fc Postnovl 
ll988l : lHils et al.lll990l ). LISA will identify several thousand 
new binaries - with the key feature of not being biased 
against short orbital periods - and will therefore contribute 
in a significant way to unveil the co mplex physical interac- 



tion of the system is not only determined by gravitational 
wave emission (and possible tidal interaction), but also by 
the mass redistribution in the system. A number of the 
currently known AM CVn syste ms are expected to serve 
as verification sources for LISA (|Stroeer fc Vecchid [iooil : 
iRoelofs et al]|2007l ). 



tions in compact object binaries (e.g. Nelemans et al.|[2004l : 
IStroeer et all l2005l : iBenacouista et all l2007l l. A particular 
class of objects are the mass-transferring AM CVn bina- 
ries (|Nelemans|[2005h . In these systems the orbital evolu- 
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In an earlier paper (jStroeer et al.ll200d ) we showed that 
in order to determine the system parameters, either the ex- 
act evolution of the gravitational wave frequency up to sec- 
ond order in a series expansions needs to be determined, 
or complementary electro-magnetic observations are needed. 
Here we revisit this result, focusing on possible short term 
variations in the equilibrium mass-transfer rate, as observed 
in many interacting binaries, and as have been proposed as 
explanation for the observed period decrease in the ultra- 
compact binary systems y407 Vul and R X J0806-I-1527 
iMarsh fc Nelemans|[2005l : iBarros et al]|2005l . and references 
therein). In addition, especially early in the e volution of AM 
CVn systems, helium novae are expected (|Bildsten et al.l 
120071 ) ■ leading to a sudden perturbation of the system. 

In Sect. [2] we describe the basic model, the way we per- 
turb the system and the calculation of the resulting gravi- 
tational wave signal. In Sec. [3] we discuss the results for one 
fiducial AM CVn system after which we conclude in Sec. U 
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2 MODEL 

2.1 Binary evolution and perturbation 

In this section we briefly review the key elements of the 
evolution of an AM CVn system, and the model adopted to 
characterize short period pertu rbations of the mas s transfer 
rate. Our formalism is based on iMarsh et al] (120041 ) to which 
we refer the reader for more details. 

The progenitor stage of an AM CVn system is char- 
acterized by the loss of orbital angular momentum due to 
the emission of gravitational waves. The lost orbital angular 
momentum drives the system to shorter orbital periods and 
smaller orbital separations. If the orbital separation is suffi- 
ciently small for the radius of the secondary star to overfill 
its Roche lobe, mass is transferred onto the primary star (ac- 
cretor) through Roche lobe overflow. The accretor responds 
by spinning up, removing angular momentum from the or- 
bit at an even higher rate. However, tidal interaction may 
force the accretor into co-rotation with the orbit, effectively 
returning the lost angular momentum. For stable systems, 
redistribution of mass in the system allows the orbit to ex- 
pand, even at decreasing angular momentum. The binary 
thus reaches a "turning point" . 

The binary under consideration is characterized by the 
masses and radii of the primary (more massive and accret- 
ing) and secondary (less massive and mass transferring) star, 
Afi,2 and i?i,2 respectively. The binary orbital separation a 
and its period P evolve due to th e change of orbital angular 
momentum Jo according to (e.g. iMarsh et al.ll20o3) 



Jo = Jgr + VGMiTJhMa + ^ (a - ^^o) 



(1) 



the three terms on the right-hand side of the equation repre- 
sent the loss of angular momentum due to gravitational wave 
radiation Jgr, mass loss at rate M2, and tidal coupling, re- 
spectively. The latter contribution is parametrized as a func- 
tion of the synchronization time-scale rs, and _Rh identifies 
the radius around the accretor with the same angular mo- 
ment um as the transferred matter (jVerbunt fc RappaportI 
1 19981 ) . fig and fio are the angular frequencies of t he accretor's 
spin and the orbit, respectively and k ~ 0.2 (|Marsh et al.l 
I2OO4I ). 

We model the mass transfer rate M2 as an adiabatic 
response of the donor 



Ma 



-/(A/i,M2,a,i?2)A' 



(2) 



w here the func t ion / (Mi, M2, a, fia) is given by Eq. (10) 



of Marsh et al.l (|2004l ) (as based on results from IWebbinM 
(|l977l )). The quantity A is the Roche lobe overflow factor, 
defined as 



A = R2- Rl: 



(3) 



where i?L is the Roche lobe's radius and a function of a 
and M2/M1 only. We use a simp l e zero -temperature mass- 
radius relation as in lMarsh et al.l (|2004l ). neglect the efl'ects 
of a realistic donor structure, that generally gives rise to 
larger donor radii, i.e. lower mass-tr ansfer rates and more 
complex mass tran sfer behaviour fe.g. lDelove &c Taamll200^ : 
iDelove et ahlbOO/l ). 

Mass transfer spins up the accretor. We therefore see 
a coupling of mass transfer rates and the response of the 
accretor due to tidal locking. This coupling can be expressed 



in a coupled differential equation system of the evolution of 
A and the evolution of fls 

1 dA Jgr kMiRf 



2R2 dt 



dQs 
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kRf J Ml 



(a - a) 



(4) 



(5) 



In the previous equation q = A/2/M1 is the mass ratio and 
the parameters (^2, Ci-l s^nd A are given, respectively, by 

d log R2 



C2 



A = 



dlogMa ' 
dlog Rl/u 
d log Ah ' 

1 + 2^^ + 



dlog k 



dlog Ml dlog Ml 



(6) 
(7) 
(8) 



We started the modelling of our fiducial AM CVn 
system in the progenitor phase and followed its evolution 
through the "turning point" towards the stable AM CVn 
sta ge assuming s t rong tidal locking (rs = 0.2yrs as favoured 
by iMarsh et al.l (|2004 ) , redbut note the discussion of the 
very large uncertainty in this value in their section 6). 
Masses for the primary and secondary star of our fiducial bi- 
nary were defined in the progenitor phase with AJi = 0.4Mq 
and M2 = O.2M0. 

We induce a perturbation by suddenly offsetting A 
(through a change in separation a that in turn changes the 
Roche lobe Rl), which directly impacts the mass transfer 
rates and spin evolution of the accreting star (see Eqs. I2I4I 
and [5|. As the mass transfer in the AM CVn systems is 
calculated implicitly from the system parameters this is the 
simplest self-consistent way of perturbing the system. The 
expected timescal e Tq for returning to equilib rium after such 
a perturbation is (|Marsh fc Nelemansll200a ) 

-I! 

We triggered the perturbation when the system reaches 
a gravitational wave frequency of 8 mHz in its stable AM 
CVn branch. We perform two kinds of perturbations. First 
we decrease the Roche lobe overflow factor, yielding reduced 
mass transfer rates, possible e.g. when nova explosions eject 
mass from the accreting star, widening the orbit. Second, 
we increase the Roche lobe overflow yielding increased mass 
transfer rates. Such an effect may be caused by suddenly ex- 
panding stellar atmospheres of the mass donor, or decrease 
in orbit if the expanding nova shell interacts with the sec- 
ondary, removing orbital angular momentum. We offset A 
such as to yield five times larger and five times smaller mass 
transfer rates (see Eq.[2]). Effectively we disturb the equilib- 
rium rates of the AM CVn system which balance Eq. |4] and 
Eq. m and force the system to evolve back to its equilibrium 
by itself. 

A factor five as perturbation for mass transfer rates is 
a quite modest change. The separation only needs to be 
changed by much less than one per cent to achieve such an 
offset. Further, we restrict this case study to a system at 
GW frequency of 8 mHz, since we refer to novae as source 
for the perturbation. The occurrence of a nova depends on 



Shortterm variations of AM CVns 3 



the ignition of the accreted layer, which happens much eas- 
ier and more often at high mass transf er rates (i.e. short 
periods) than later (jBildsten et al.l[2007l ). 



2.2 Gravitational wave signal 

We model the GW that is emitted from such a system at 
the leading Newtonian quadrupole order, only keeping the 
+ polarisation and working in the real domain. This yields 
for the GW strain in the solar reference frame 

hawit) = A{t) cos(cE'Gw(t)) = A{t) cos{2nt few (t)), (10) 



with 



1 2Mi{t)M2{t) 



D c4 



a{t) 



(11) 



the amplitude of the strain (Afi,2(t) masses of the stars in a 
binary as affected by mass transfer, D, a{t) distance to and 
varying separation of the system respectively, G, c Keplers 
constant and the velocity of light respectively, $Gw(i) the 
phase and few the frequency of the gravitational wave). We 
calculate the period of the AM CVn according to 



(P(t)/27r)' = a{tf/G{Mi{t)+M2it)) 



(12) 



(Kepler's third law) to derive the gravitational wave fre- 
quency 



fGwit) = 2/P{t) 



(13) 



(leading order quadrupolar radiation as found in circular 
orbits, we drop the subscript GW hereafter). 

We do neglect the LISA instrumental response function 
in this research, and thus remain in the solar reference frame 
through out. We do note that this does not have repercus- 
sions for the conclusions, as the response function does not 
change the GW frequency evolution (it merely modulates 
it). In the remainder of this paper we drop the explicit ref- 
erence to the time dependence in the parameters. 

We expanded the evolving phase of gravitational wave 
emission from the AM CVn system as a Taylor series up to 
second order to determine the influence of the perturbation 
to the GW signal. 



(14) 



<i>Gw{t) = 27r/ot + nfl,'h^ + ^f^^h^ 

We defined with /o, /o^' , /o^' spin-down parameters for the 
system, here stating the gravitational wave frequency, its 
first and second derivative in time at the origin of the series, 
time nominally set to zero at an instant after the sudden 
offset of A. 3>Gw(i) is further the number of wave-cycles 
one expects to record at LISA within a given observation 
time. For this the running time t is replaced with the total 
observation time T. In this view, the three terms to ^gw{T) 
are contributing separately to the total number of wave- 
cycles. 



3 RESULTS 

In Fig.[T]we show the frequency evolution of the unperturbed 
system (solid line) and the perturbed system (dashed and 
dotted lines) . We observed a change in drift of the frequency 
of gravitational wave emissions after a perturbation struck 



7.995 - 



7.9S5 



7.ga 



smaller mass transfer 



larger mass^-. 
trarvsfer 



100 



200 300 
t [yr] 



400 



500 



Figure 1. The evolution of the gravitational wave frequency of 
the fiducial AM CVn system over time, with the x axis counting 
from the beginning of the perturbation. The solid line denotes the 
equilibrium evolution, the dotted line an evolution with a mass 
transfer rate five times smaller than in equilibrium at onset of 
perturbation, the dashed line an evolution with a mass transfer 
rate five times bigger than in equilibrium at onset of perturbation. 



the system. Perturbations that lead to increased mass trans- 
fer rates forced the system to separate at a higher pace, leav- 
ing the gravitational wave frequency to decrease at a faster 
pace as compared to a secular evolving AM CVn at the same 
stage of evolution (compare Fig.[T]dashed line with the solid 
line). Perturbations that decreased the mass transfer rate, 
in contrast, lead to a short in-spiral of the system over 59 
years (see Fig. [T]dotted line). The system experiences a new 
"turning point" and then evolves to lower frequencies again. 

We determined the total number of wave-cycles as de- 
tected during a one and five year mission of LISA and the 
contributions from the spin-down parameters to this total 
number by "starting" an observation at each individual year 
of the shown frequency evolution of the AM CVn in Fig. [T] 
(see Fig. [2] and Fig. |3J. We found, that the shapes of the 
curves in Fig. [2] and Fig. |3] do not change between a one 
year observation and a five year observation, they simply 
scale according to and T"^ for the first and second deriva- 
tive respectively. We therefore show explicit examples for a 
5 year mission only. 

The number of wave cycles contributed by /q^' clearly 
mirrors the altered frequency evolution. Perturbations that 
force increased mass transfer rates (dashed lines) yield faster 
separating systems, i.e. large negative contributions by /q^'. 
We find the number of wave cycles from the first spin-down 
parameter to increase from -123 wave cycles to -990 wave 
cycles when integrating the number of wave cycles over a 
five year observation starting an instant after the sudden 
offset of A. Perturbations inducing decreased mass transfer 
rates (dotted lines) yield a short in-spiral of the system, we 
thus see first moderate positive then moderate negative con- 
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Figure 2. Number of wave-cycles as contributed by the first 
derivative of the gravitational wave frequency (spin-down param- 
eter) over time as seen by LISA during a five year observation. 
Time as shown on the x axis refers to the timeline since pertur- 
bation struck the system, and a hypothetical LISA observation 
is started at every instance of this timeline to yield shown wave- 
cycle evolution. The solid line denotes the equilibrium evolution, 
the dotted line an evolution with a mass transfer rate five times 
smaller than in equilibrium at onset of perturbation, the dashed 
line an evolution with a mass transfer rate five times bigger than 
in equilibrium at onset of perturbation. The inlay shows a zoom 
in in order to cover the full spread of the evolution with increased 
mass transfer rate at onset. For details see text. 

tributions by fg^^ as separated by the new turning point in 
frequency. Here the significant alteration is more the reverse 
of the drift of the system, with -1-50 wave cycles immediately 
after the impulse instead of the quoted -123 wave cycles. 

(2) 

Interesting results can be found for /g . This spin-down 
parameter normally does not contribute to the total num- 
ber of wave-cycles for a system at 8mHz even in a five year 
observation. In case of a faster separating system however 
we suddenly see a significant increase in the contribution to 
the total number of wave-cycles, in fact by a factor of 10*^. 
Suddenly the second derivative becomes visible within a five 
year mission with +76 additional contributed wave-cycles, 
compared to a null contribution beforehand. In case of the 
system with decreased mass transfer we see an increase in 
the number of contributed wave cycles, but still the contri- 
bution is too small to be detected. 



4 DISCUSSION AND CONCLUSIONS 

We studied the reaction of the frequency evolution of a fidu- 
cial AM CVn system with a GW frequency of 8mHz, to 
sudden, but modest, perturbations of its orbital parameters. 
These perturbations are inspired by the observed and theo- 
retically expected short term variations due to e.g. nova ex- 
plosions. The system remains bound, changes its mass trans- 
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Figure 3. Number of wave-cycles as contributed by the second 
derivative of the gravitational wave frequency (spin-down param- 
eter) over time as seen by LISA during a five year observation. 
Conventions as in Fig. |2] 

fer rate, but on a time scale of hundreds of years, evolves 
back to its equilibrium mass transfer rate, as expected. How- 
ever, if any of the thousands of individual binaries that LISA 
is expected to follow during its life time, is in such a per- 
turbed state, the measured frequency evolution can be dra- 
matically different than for the equilibrium case. The sign 
of the frequency derivative may change and the number of 
cycles contributed by the first and second derivative may 
increase by several tens to several hundreds in a 1 or 5 year 
mission life time. For template data analysis methods, the 
detection of these systems should still be secure, if sufficient 
spin-down parameters are used (up to second derivative, es- 
pecially for a 5 year mission). However, care should be taken 
with deriving system parameters purely on the GW data, as 
the perturbations can give rise to wildly different derivatives, 
which, if interpreted as equilibrium evolution, would lead to 
extreme mass estimates, or for the systems with reversed 
evolution to t he conclusion that no mass transfer is present. 
In addition in IStroeer et al] (|2005h we concluded that /q^' 
is needed for distance estimations, but this works only when 
the system evolves in equilibrium. 

It is difficult to estimate the number of LISA bina- 
ries that are affected by such short term variations. Many 
observed interacting binaries seem to have short term pe- 
riod evolution differ ent from the secular evolution (e.g. 
iBaptista etall (|2003l )). suggesting it may be a quite com- 
mon phenomenon. For the helium novae in AM CVn systems 
we can make a simple estimate. At the shortest periods, i.e. 
high mass transfer rates and for high accretor masses the 
amount of accum ulated mass needed to ignite the nova is 
about 0.001 Msun (|Bildsten et al.ll2007l ). For a mass transfer 
rate of 10~* Msun/yr, this leads to a nova every 100 years, 
influencing the evolution essentially all the time so we expect 
several LISA sources to be affected. However, most of the 
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LISA sources will have evolved to lower mass transfer rates, 
where the nova frequency is much lower and the chances of 
catching a system in the recovery phase are low. 

We also tested repercussions of short term perturba- 
tions at an evolutionary stage of the fiducial AM CVn at 
which the orbital separation yields GW frequencies of order 
3 mHz. This stage is set 3.5 million years after we found our 
AM CVn at 8mHz. Overall, tidal coupling is much less at 
such a separation so that the response to a shortterm per- 
turbation is less dynamic. We find the timescale of the per- 
turbation, Tc, to be 691 yrs compared to 47.9 yrs. Therefore 
the maximum amount of increase in wavecycles, here found 
in the increased mass transferring scenario, is only a mere 
2.5 wavecycles, and the second derivative is never visible in 
case of increased and decreased mass transfer perturbations 
at proposed scale. However, decreased mass transfer is still 
able to trigger a new turning point in the frequency evolu- 
tion of the binary, the system inspirals again over a period 
of roughly 2000 years. We therefore conclude that although 
the magnitude of wavecycle increase is minor at 3 mHz the 
actual ability of the system to shortly inspiral is a reason of 
concern once again for parameter inference. 

W e thus re-enforce our earlier conclusion (|Stroeer et al.l 
I2OO5I ) t hat complement ary electro-magnetic data from e.g. 
GAIA (|Nelemansll2"009l ) is an important ingredient for the 
disentangling of different processes that may contribute to 
the frequency evolution of compact binaries, such as mass 
transfer, tidal interaction and short term variability and the 
determination of accurate system parameters. 



B. E., Steeghs, D., Morales-Rueda, L., Marsh, T. R., & 
Nelemans, G. 2007, ApJ, 666, 1174 

Stroeer A., Vecchio A., Nelemans G., 2005, Astrophys. J. 
Lett., 633, L33 

Stroeer A., Vecchio A., 2006, Classical and Quantum Grav- 
ity, 23, 809 

Verbunt P., Rappaport S., 1998, ApJ 

Webbink R. F., 1977, Astrophys. J., 211, 486 



REFERENCES 

Baptista R., Borges B. W., Bond H. E., Jablonski P., 
Steiner J. E., Grauer A. D., 2003, MNRAS, 345, 889 

Barros S. C. C, Marsh T. R., Groot P., Nelemans G., Ram- 
say G., Roelofs G., Steeghs D., Wilms J., 2005, MNRAS, 
357, 1306 

Benacquista M. J., Larson S. L., Taylor B. E., 2007, Clas- 
sical and Quantum Gravity, 24, 513 
Bender P. L., 1998, Bulletin of the American Astronomical 

Society, 30, 1326 
Bildsten L., Shen K. J., Weinberg N. N., Nelemans G., 

2007, ApJL, 662, L95 
Deloye, C. J., Taam, R. E., Winisdoerffer, C, & Chabrier, 

G. 2007, MNRAS, 381, 525 
Deloye, C. J., & Taam, R. E. 2006, ApJL, 649, L99 
Evans C. R., Iben L J., Smarr L., 1987, ApJ, 323, 129 
Hils D., Bender P. L., Webbink R. P., 1990, ApJ, 360, 75 
Lipunov V. M., Postnov K. A., 1988, APSS, 145, 1 
Marsh T. R., Nelemans G., 2005, Mon. Not. R. Astron. 

Soc, pp 837— H 
Marsh T. R., Nelemans G., Steeghs D., 2004, Mon. Not. R. 

Astron. Soc, 350, 113 
Nelemans G., 2005, in ASP Conf. Ser. 330: The Astro- 
physics of Cataclysmic Variables and Related Objects AM 
CVn stars, pp 27—1- 
Nelemans G., 2009, Classical and Quantum Gravity, 26, 
094030 

Nelemans G., Yungelson L. R., Portegies Zwart S. P., 2004, 

MNRAS, 349, 181 
Roelofs, G. H. A., Groot, P. J., Benedict, G. P., McArthur, 



